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ABSTRACT 


A form of the numerical lunar theory is suggested for plane- 
tary satellites moving in orbital planes with higher inclinations 
relative to the orbit plane of the Sun. The number of terms given 
explicitly in final solutions is too small to warrant the noncritical 
application of the written portions of the series to the cases of 
higher inclinations. In addition, the convergence of these series 
is slow near resonance angles. Convergence problems can be 
remedied in part by arranging the theory in such a way that the 
squares of the mean motions of the perigee and node are deter- 
mined instead of the first powers of these motions. BrendeTs 
coordinates are used, and the equations of motion are reduced to 
the problem of two disturbed linear oscillators. Following 
Hansen the numerical values of all the basic parameters of the 
lunar theory are substituted from the outset and the process of 
iteration is used to solve the problem. The mean orbital plane 
and the constants of integration are chosen in such a way that the 
expansions of the coordinates contain only purely periodic terms. 
The polar angle v of the projection of the satellite on the mean 
orbital plane is the independent variable of the theory. With this 
choice of the independent variable, the differential equations of 
the problem become the equations for linear oscillators. 
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A NUMERICAL THEORY OF SATELLITES IN BRENDEL’S COORDINATES 


by 

Peter Musen 

Goddard Space Flight Center 


INTRODUCTION 

In this paper, we suggest a form of the numerical lunar theory that can be used for planetary 
satellites in orbits with high inclinations to the orbit plane of the Sun relative to the planet. The 
classical analytical lunar theories by Delaunay (1860), Hill (1878), and Brown (1895) assume the 
smallness of the inclination I, and sin 1/2 is considered to be one of the small literal parameters 
with respect to which the coordinates or the elements are being expanded. 

Despite the elegance of these literal developments, the number of terms given explicitly is too 
limited to warrant an uncritical application of the written portions of the series to the cases of 
high inclinations. The theory of Delaunay— and, one might suspect, the theory of Hill and Brown- 
have resonances near I ~ 39°, 141°; this affects their applicability for higher inclinations, because 
of the slow convergence of the expansions near these critical inclinations. This slow convergence 
can be partly remedied by determining the square of the mean motion of the pericenter instead of 
the motion of the perigee itself, because the square of the mean motion of the pericenter can be 
written in a more closed form than can the motion of the pericenter itself. Furthermore, the prob- 
lem of convergence can be removed to a considerable extent if numerical values of the parameters 
are substituted from the outset — especially the parameter m, the ratio of mean motions, or its 
equivalent. This fact is well known from the Hill- Brown lunar theory, in which the numerical value 
of m is used. We suggest here the numerical substitution of all small parameters and a purely nu- 
merical theory. In this respect we decided to follow in Hansen's footsteps (1862). 

An important advantage of the purely numerical approach is that the terms are computed in 
large groups at each iteration step. The number of cycles in the iterative process depends upon 
the values of the parameters and upon the final accuracy needed for computing the ephemeris of 
the satellite. In every lunar theory the small divisors and the long-period effects constitute a 
difficult computational problem. The only safeguard against distortion of the critical terms in a 
numerical theory is computation with a formal accuracy higher than that needed in the coefficients. 
In other words, one must program with tT double precision.” 

We suggest also the use of BrendeTs coordinates s, £, and </>, where s is the perturbed value 
of the projection of the unit vector of the satellite on the Laplacian vector, £ is the elevation of the 
satellite relative to its mean orbital plane, and ^ is proportional to the areal velocity of the projec- 
tion of the satellite on its mean orbital plane. The mean orbital plane and the constants of integra- 
tion are so chosen that the expansions of BrendeTs coordinates contain only purely periodic terms. 
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Then the expansions of the satellite coordinates in the inertial system of coordinates will also 
consist of purely periodic terms. 

Some points of affinity between Hansen's theory and the theory in Brendel's coordinates should 
be pointed out. Hansen's lunar theory makes use of the osculating orbit plane as a reference plane. 
The equations of motion of the satellite in the reference frame rigidly connected with the 
osculating orbit plane have the same form as in the inertial frame. For this reason a frame of ref- 
erence rigidly connected with the osculating orbit plane was termed ideal by Hansen. Its intro- 
duction facilitates treatment of the perturbations of the satellite in the osculating orbit plane but 
requires determination of the position of the osculating orbit plane relative to the inertial system 
as an intermediary step before the expansion of the coordinates is done. In our theory we choose a 
reference frame with the<f- and 77 -axes in the mean orbital plane and with the z-axis normal to this 
plane in such a way that the equations of motion of the satellite relative to the frame (£ 77 z ) will differ 
in form by the terms of higher orders from the corresponding equations of motion in an ideal sys- 
tem. As compared with Hansen's theory, such a choice of the reference frame, together with the 
use of Brendel's coordinates, permits us to obtain the rectangular coordinates in the inertial 
system in a more direct way, bypassing the determination of the position of the osculating 
orbit plane. 

We choose the polar angle v of the projection of the satellite in the (£77) system to be the in- 
dependent variable of our theory. With this choice the differential equations for s and £ take the 
form of the differential equations of disturbed linear oscillators. These equations resemble the 
corresponding equations of Brendel's planetary theory (1925) but contain also some additional 
terms which reflect the motion of the "nearly ideal" system (^z ) relative to the inertial system. 

The form of the differential equations of the lunar theory given here also favors the applica- 
tion of the Krylov-Bogolubov method (1961), if a purely analytical form of the solution is preferred. 
The expansion of Brendel's coordinates is obtainable in the form of trigonometric series in four 
arguments Wj , w 2 , w 3 , and w 4 , each linear in v. The angle is the mean argument of the latitude, 
w 2 is the mean true anomaly, w 3 is the mean (relative to v ) longitude of the Sun, and w 4 is the 
mean longitude of the ascending node. 

The form of the differential equations for s and £, resembling the differential equations for 
oscillators, suggests the determination of squares of the frequencies of w x and w 2 rather than the 
frequencies themselves. The positive side of this is that the squares of frequencies have a more 
closed form relative to the mean inclination than do the frequencies themselves. We also might 
add that Brendel's planetary theory, as well as the lunar theory in Brendel's coordinates, re- 
sembles the Laplacian lunar theory (with the requirement for smallness of the inclination removed). 
It seems that in our time we are witnessing a revival of Laplacian ideas in more or less hidden 
forms. 

There is a tendency now to represent the coordinates or the osculating elements in terms of 
arguments other than time. Thus, in the case of an artificial satellite, the perturbations of the 
elements are obtained in terms of the true anomaly and also in terms of the arguments which are 
linear with respect to time (Kozai, 1959, and Brouwer, 1959). In Musen's (1961) theory, the true 
orbital longitude is taken as the independent variable. There is a similar tendency to represent the 
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perturbations of the elements of the satellite in terms of the true, or of the eccentric, anomaly of 
the satellite and in terms of the true anomaly of the disturbing body. All these trends belong to 
the circle of Laplacian ideas. 

BASIC DIFFERENTIAL EQUATIONS 

Let us assume that the trajectory of the Sun relative to the planet is a fixed ellipse. This also 
means that the satellite mass is considered negligible. We put the X- and Y-axes into the plane of 
the solar ellipse; the Z-axis will be normal to this plane. The system (XYZ) is the inertial system 
to which the motion of the satellite is referred. 

The satellite’s mean orbit plane (xy) in the present theory has a constant inclination I relative 
to the XY-plane and rotates around the z-axis at such a speed that it approximates the osculating 
orbit plane with an accuracy up to the periodic perturbations in the longitude of the ascending node 
and in the inclination. Furthermore, all frequencies of the arguments of the theory, as well as the 
constants of integration, must be so defined that no purely secular or mixed terms appear in the 
expansions of the satellites coordinates relative to the system (xyz). Let 

y ~ yk 

be the angular speed of rotation of the xy -plane, Pbe the unit vector along the line of the ascend- 
ing node of xy on the XY-plane, R be the unit vector normal to xy in the direction of the angular 
momentum, and 


Q = R x p . 

We place the axes x and y along P and Q, respectively; the z-axis is along R. We have 

y = \Q + v R , (1) 


where 


K - y sin I , 


y cos I 


( 2 ) 


The differential equation of motion of the satellite relative to the system (xyz ) has the form 


d 2 7 _ dr dy ___ r 

+ 2y X -fa + X r + y x (y x r) - - ~ + grad- Q 


( 3 ) 


We decompose the position vector r into the component p 9 which is the projection of r on the 
xy-plane, and the component z 5 normal to that plane; we have 


r - p + z R 


( 4 ) 
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Substituting Equations 1 and 4 into Equation 3, gives 


/ dp dv \ / dz d\ \ p 

x dt + dt Pj + aF + dF j P = " f: 

d 2 z _♦ / dp d\ _A z dF 

dt 2 P ’ ^ dt + dt Pj r 3 + dz 


+ grad-F 


where 


F = j [y 2 x 2 + (vy-\z) 2 ] + H . 

The scalar differential equations of motion as obtained from Equations 5 and 6 are 


■ dy 

2^ dF + 
< 

dv 

dt 

yj + (f k Ht 

dk \ 

+ dF A 

= - 

d 2 y 

( 

dx dv \ 

y 

dF 

dt 2 + 

\2v 

dt + dt x ) 


+ dy 

d 2 z 


dx dk \ 

z 

dF 

dt 2 

dt + dt x ) 


+ dz 


2L 

r 3 + dx 


Equations 8 through 10 are almost identical in form to those suggested by Brouwer (1958) for the 
development of the Hill-Brown theory of motion of artificial satellites. In Brouwer f s work k and 
v are constants, and thus the terms involving dk/dt and dv/dt are absent. 

Equations 5 and 6 can be transformed to a form closely resembling the equations of motion in 
an ideal system and more convenient for the determination of Brendel's coordinates. We refer the 
motion of the satellite to the system («fpz) rotating relative to (xyz) with angular velocity -^R. In 
order to transform the equations of motion from the system (xyz) to the "nearly ideal TT system 
(£r)z ) we must replace, in Equations 5 and 6, 


dp , dp 

hf b y dF - v R * p 


- dp dv -* _ , 

2 vRi< dt"dt E,< ^' 1 ' v 
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in accordance with the basic kinematic rules for the relative velocity and relative acceleration. 
The derivatives 


dp d 2 p 

dt and dt 2 

now mean the velocity and acceleration relative to the system (£, 77 , z). We obtain 

d 2 ~p 
dt 2 


d 2 z 
dt 2 

where we put 

$ = ^ k 2 (x 2 + z 2 ) - \vyz + fi . (13) 

The form of Equations 11 and 12 now resembles, as closely as possible under the circumstances, 
the form of the equations of motion in the ideal system of coordinates. 

In the formation of the differential equations for BrendeTs coordinates, it is convenient to 
choose the angle v between p and the £-axis as the independent variable of our theory. This seems 
to be a most natural choice, because v represents the closest analog of the true orbital longitude. 
In accordance with this choice, we set the angular speed of rotation of the (xy)-plane to be uniform 
with respect to v. In other words, we set 


~p ( dz dA. \ — 

+ grad^O - (2^ dt + dt Z J P ’ 

/ dx d\ \ 

^ dt + dF X J • 


z dO 

+ “ 5 — + 

r 3 dz 


( 11 ) 


( 12 ) 


V 


dv 

(1-8) 3t 


where g is a constant. Then we have 

dv 

K ~ (i - g) dt tan 1 ■ 

dv 

y - (1 - g) sec I dt * 

The angular distance w of p from P is obtained from the equation 

w = V - fi'dt - gv - & 0 , 
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where 6 0 is a constant. This angle is one of the basic arguments in the expansions of Brendel T s 
coordinates into trigonometric series. 

^ From Equation 5 we obtain, by projecting the left and the right sides on the directions of p and 
of R xp: 


d 2 p /dv\ _ p <9$ / dz dk\ 

~fa~2 ~ P \dt) " r 3 + dp “ l 2 ^ dt + z dt J 

d / dv\ _ <9$ / dz dk\ 

dF \P 2 di) ~ <9w + (2^ dF + z dF/ <° sin ' 

and, in addition, from Equation 6: 

d 2 z z (90 (dp cos w dA. \ 

"^7 = ” + aF + ( di~ + dt p cos H 


(14) 

(15) 


(16) 


At this point we introduce, instead of p, v, z, Brendel coordinates *A, s, £ by means of the equations 


2 dv 
p dt 

= Vp . 

(17) 

P 0 ( x " 0) • 

p 0 - constant , 

(18) 

P 0 C 1 - 

p 


P ~ 1 +s 

- Its ’ 

(19) 

z r 

■ 

(20) 


It must be emphasized that p is not the osculating semi-latus rectum in the sense of the theory of 
variation of astronomical constants, because the system (<fr?z) is not ideal. However, we can put 
p = p 0 (1 where p 0 is a constant and 0 is the perturbations whose expansion contains only 
purely periodic terms. In Equation 19 s is analogous, but not identical, to the osculating value of e cos f. 

Pursuing this development further, we take v either in the form 


V 


( 1 - 



( 21 ) 


or in the form 


V 


p<r 3 / 2 (i-g)a-'/')' 3 / 2 (i +s > 2 . 


( 22 ) 



and we make use of the relation 


K - v tan I . 

From Equations 15, 17, 20, 21, and 23 we deduce 


d /p 

dv 


p 2 <90 

yp <9w + (* 6) 


ip 


d log p 
dv 


*»&*< 


d log v\ 

dv J 


sin w tan I 


Making use of Equations 18, 19 and 22 we obtain from the last equation, 


d yp 

dv 


2 p z 

P 0 <9w + 


(1 


d0 dM 

■ g) £ d^ “ 4 C1 "'/O dTl 


s in w tan I 


(23) 


(24) 


From Equations 19 and 17 we have 


dp _ 1 A + s d ip d s \ 

dt " ~ ^ yi - 0 dv + dvy 


(25) 


Differentiating Equation 25 and taking Equation 17 into consideration again, we obtain 


p 


2 


d 2 P 
dt 2 


d 2 s 1 -t s d 2 \p 1 ds dip 3 1 + s /di/\ 2 

“d^ 2 ” dT 2 " "2 1'Hv dv _ 2 (1 ^ )2 V^7 


(26) 


The two last equations appear in BrendeTs work. However, the meaning of v in Brendel T s equa- 
tions is different from the meaning assigned to it here. 

From Equations 14 and 17 we have 


d 2 p p /p \ 3 „ <90 /_ dz d log u\ _ 

~ p " + P ~d^ ~ v ( 2 dv + z dv ) yp cos w tan I 


(27) 


and taking Equations 19 and 21 into account we can rewrite this equation in the form 


d 2 p p /p \ 3 0 <90 ( d£ 

~^2 = P- [T) + P 2 ^-(l-g)(l+s) 


2C 


d log p d log 


dv 


+ i 


dv 




cos w tan I 


Substituting Equations 19, 22 and 


r - p^ 1 + £ 2 


(28) 
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into the last equation, we obtain 


fj j_P = 

dt 2 


s + 1 - (1 + £ 2 )~ 3/2 + p 2 |“- (l-g)(l +S) ^2 dv - 2 


cos w tan I 


(29) 


Comparing Equations 26 and 29 we deduce 


d 2 s 2 <9$ 

— + s = A + B - p -Qp 


(30) 


where 


1 + s d 2 0 3 1 ds d0 

1 ~ 0 ^ v 2 2 1 “ \p dv dv 


(31) 




di I £ dA _ 

• - - - . . l cos w tan I 


(1 -g) (1 + s) 2 d7 dv 


(32) 


Next, from Equations 16, 17, and 21, 


d 2 £ 2_dp d£_ £ d 2 p 

dt 2 + p dt dt + p dt 2 


_L 1 dO 

r 3 + p dz 


P / d log p d . , 

+ — (l“g) (2 cos w 2g sin w + cos w — — 1 tan I. 




Substituting Equations 19 and 22 into the last equation gives 


d 2 j 2dpd$ id = _ 150 

dt 2 + p dt dt + p dt 2 r 3 + p 5z 


p (\ 1 d0 

" (1 ' g) ^\2 T ~ r ^' d 7 cosw 


+ 


2g sin w] 


tan I 


(33) 


Differentiating 


d£ = j^dS 
dt p 2 dv 


(34) 
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gives, after some easy transformations, 


d 2 g _ p_ fd 2 C 1 d£ d log p dZ d logp\ 

dt 2 p 4 \dv 2 2 dv dv “ 2 dv dv ) 


(35) 


Substituting Equations 18 and 19 into this equation gives 


d 2 j _ _p_ fd 2 £ _3 1 d0 d£ 2 ds dA 

dt 2 p 4 \dv 2 + 21-0dvdv + l+sdv dv^ 


(36) 


Multiplying Equations 25 and 34 together gives 


_l dp d£ _ ( l d 0 d£ 1 ds d£ 

p dt dt ^4 yl ~ 0 dv dv + 1 + s dv dv 


(37) 


From Equation 27, taking Equation 21 into account, 


ldV 

1 P 


( dz 

d log u\ 

cos w tan I 

P dt 2 

r 3 p 4 

P dp) 

- (i - g) (2 a v + 

z dv J 

P 


or 


A d2 P 
p dt 2 



(1-g) 


dt 

dv 


+ u 


d log p d log v 

dv + ^ dv 


cos w tan I 


and finally, eliminating £ and v in favor of s and 0 as before, 


_1 d2 P 
P dt 2 



P 3 ao\ 

p dp J 


' d-g) 


L I 

\^ 2 dv 2 


_L_ d ^| 

1-0 dvy 


cos w tan I 


(38) 


Substituting Equations 35, 37, and 38 into Equation 33 gives 


d 2 £ 

+ l +2g(l - g) sinw 


P + 


£l (i®. 

P y?z 



(39) 


where 


„ _ 1 1 Ap d£ fdl 2 1 1 + l 2 d A 

^ 2 1 - 0 dv dv + ( 1 8 ) l dv 2 1 - 0 dv J cos w * an ^ 


( 40 ) 
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If in the development of $ we eliminate z in favor of £ by means of Equation 20, then 


d$ _ d$ 

d z and dp 


in Equations 30 and 39 will be replaced by 

1 d<t> d<P 

p <9£ and dp ' p <9£ * 

respectively, and the complete set of differential equations of the problem is 


d 2 s 

dv 2 


+ 


s 


A + B - p 


f <90 
PJp “ 



d 2 £ 
dv 2 


+ £ + 2g( 1 - g) sin w 



dv 


2p 2 

p 0 


+ d-g) 


s dv 


- 4( 1 - i/>) 


di 

dv 


sin w tan I 


dt 

dv 


= p 0 3/2 (l-^) 3/2 (l + s)- 


where A, B, and P are defined by Equations 31, 32, and 40, respectively. 
Combining the equations 


dv ; _ , 

at _ yc 1 +m ) p 0 . 


1 +s' 


s e cos v 


with Equations 17 through 19 gives 


dv' 

dv 


= m 0 - <j/) 


3/2 


(W 


where 


- — (T 


(30') 


(39’) 


(24') 

(41) 


(42) 
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We now split the right side of Equation 42 into a constant part m, which is of the order of m 0 , and a 
purely periodic part Q. We thus have 


v' - w 3 + Sv' 


(42’) 


where 


w 3 - mv - v Q ' , v Q ' - Const 


(43) 


and 


Sv' 



(44) 


The integral in Equation 44 is taken in a formal way, without adding any constant of integration. 
The angle w 4 is one of the basic arguments of the theory, analogous to the argument mv + const of 
the Laplacian (1802) theory. However, the meanings of v in the theory presented here and in the 
Laplacian theory are different. The remaining arguments are 

W = CV - 7T 0 , 

which is introduced in the process of integration of Equation 30% and 

W 4 = gv - 0 0 , 


the mean longitude of the ascending node of the satellite. We shall see in one of the next sections 
how these four arguments appear naturally in the expansion of the disturbing function. 


THE GENERALIZED CAUCHY NUMBERS 

The Cauchy (1885) numbers N_ p . can be defined as the contour integrals 

= ” " ( z **)’ ( z * *)’ dI 


(45) 


taken along the unit circle | z| = 1. Cauchy has shown that their introduction permits one to ob- 
tain in a simple way some basic developments associated with the elliptic motion. 
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In the development of the higher harmonics in the solar disturbing function into trigonometric 
series for a satellite in a highly inclined osculating orbital plane, we found it convenient to intro- 
duce the set of rational numbers 


N 


V 

-p. i .q 


1 

27Ti 



z _p_1 cy 




(46) 


whose properties resemble closely the properties of the Cauchy numbers. The symbol C. v designate 
a Gegenbauer polynomial; p, j , q are integers; j and q are nonnegative; and p can be positive, 
negative, or zero. The generalized Cauchy numbers NZ pJ represent a useful instrument in ob- 
taining the disturbing function of the lunar theory by means of electronic computers. 

Writing the classical expansions of Gegenbauer polynomials (Erdely et al., 1953) in the form 


C y 


hi) 


1 y-> (V, m)(v, j -m) / 1 \ 

2 + zi -2 m ) 


z + M Ciza] 


c v 


' (-!)■" (v, j ~m) 1 / J.V" 

/ . ">! ( j - m)! 2 >~ 2m \ z ) 


r v r + z \ - . (V, j) y (-j. m) (j + V, m) 1 / 1\ 2 

” \ 2 ) ~ ' ' ^ L n ' * V 

\ / m=o 1 m I ( l t m) 


C* 


2 j \ 2 


(2v, 2j) (-j, m) (j + v, m) (~l) m 

2 


y (-j. m] 

(2j)! Z- / 1 ) 

m- o ( v + * 2 , m ) ( 1 , 


m) 




/z + 


CZ 


2j + i \ 2 


. - y 0> j + 1) V 1 (-j . m) ( j + V + 1, m) 1 / l\ 2m+1 

(-1) j ! Z. /3 \ + V) 

m= o I 2 ", m (1, m) 


f z + 


CZ 


2j + i \ 2 


(2j + 1)! 


(2^’ 2j + 1) y-> (~j, m) (j + 1/ + 1, m) (-1) 
I—i I 1 \ 

m = o I v + "j * ml (1 , m ) 


2m+l 


H)H)' 


and substituting them in Equation 46, we obtain a set of formulas giving the numbers NZ p . in 
terms of the standard Cauchy numbers: 


N v = if 1 (?* m ) (^» J ~ m ) / v 

"P.j.q 2 / m! ( j " m)! \ W “P“j + 2m l 0, q + W -p4j-2m, 0, qj ’ 


(47) 
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( 48 ) 


(i/d 


N v = y 1 

-p. j. q / rn! (j -m)! " 2 j-2m -p, j-2m, 


N*' 


p ,2 j.q 


(-iy 


(M ) m)(j +v. m) 1 


02 y (-j. 

M L(± 


A \ „ 

o C 1 ’ 


m) 


22 m -p, 2m, q ’ 


N v 


~P.2j.q (2j) 


- (2^. 2j) y(-j, m)(j + v, m) (~l) w 


■) ) y (-] , m, 

r LLi 


Fo 1 (li m ) 


22 m -p, 0, q+2m ’ 


= /_,sj ( v » i + 1 ) y (-J. rc)(j + v + l, m) 1 


N v = ...j o. 1* 1 ,) y 

N -p.2 i + l,q ( 1) j! ^ 


N v 


(| (1, m) 

(2v, 2j + 1) V~i (~j, m)(j + + m) (-l) m 


22 m -p, 2m+ 1, q 


(-j . rn J 

p,2j+l,q ( 2 j +1)! / / 1 \ 

m=° tv + 2 ", ml (1, m) 


22 m+l ~p, 1, 2m+q 


From these formulas it is evident that, as is the case for Cauchy numbers, 


N v . = 0 

~P. J, q 


if j + q - p is negative or odd and 


N u 

p. . 


= (~l) q N v . 
q v ' ~p. i, 


(49) 


(50) 


(51) 

(52) 


(53) 


(54) 


In particular, 


(55) 

(56) 


A set of recurrence formulas can be established between the generalized Cauchy numbers. 
Making use of the recurrence formulas 

C i +2) c/; 2 (X) = 20 + j + 1) X c; t] (X) - (2v + j ) C? (x) , 


jc; (x) = 2v [x Cjl* 1 (x) - C,17 (x)] 
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and setting 


(j +2v)C/' (x) = 2v [C/'* 1 (x) - xCjlV (x)] , 

i C. v (x) = ( j - 1 + 2v) x C*. t (x) - 2v(l - x 2 ) C (x) , 



we deduce the recurrence relations 


N 


V 

-p, j + 2, q 


V + j + 1 
j + 2 


(n 


V 

~p + Ir 


i + l. q 


N - V P - 


l. j + l. q ] 


2v + j 
' J + 2 


N. 1 ' 


(57) 


N 


V 

"P- 


j. q 


v 

T 


(N v + 1 
V “P+ 1. 


+ N v+l 

J - 1. q “ P" 1. I “ 1. q 


- 2N 


v + 1 

-P, j “ 2, q 


) 


(58) 


j + 2v 




N V+1 

-P+ 1, j - 1. q 


N v + 1 . , ) 

■p i. j ■ l. q^ 


(59) 


N v 


1 + 2 ^ 
2j 


(n v + N v \ 

V -p+1. j-1. q 1 -P-1, j-i. qj 


V 



N 


v- 1 

“P d-2 ,q + 2 


(60) 


which, combined with Equations 53 and 54, can be used to compute the table of the generalized 
Cauchy numbers. 


A TRANSFORMATION OF THE ADDITION THEOREM FOR GEGENBAUER POLYNOMIALS 

We shall transform the addition theorem for Gegenbauer polynomials to a form suitable for 
obtaining the expansion of the disturbing function of satellites in orbits with high inclinations. 
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The transformed addition theorem shows that the generalized Cauchy numbers represent a useful 
instrument for obtaining the expansion of the disturbing function, especially of the terms of higher 
order. We have 

n 

(cosH) = ) 2 2m k v sin m 6 C l/ + m (cos 6) sin m \pC ! *l +m (cos <//) C v ~ 1/2 (cos c£) , (61) 

q x 7 / n t m n ~ m v ' n m N m N ' ' ' 

m= 0 

where 


cos H - cos 9 cos *// + sin 6 sin Kp cos 0 


(62) 


and 


k v 

n , m 


(n-m)! 


r(2 v-i) ■ r 2 p+m) 
T 2 (v) •P( 2 v+n +m) 


(2 v + 2 m - 1 ) . 


(63) 


We set 


x - exp \Q , y exp it// 


(64) 


Then we have the expansion of the form 


j , k~+n 


C n v (cos H) - [n, j , k] v x j y k 


j , k = ~n 


(65) 


where [n, j, k]^ are polynomials in cos 4>. We have from Equation 65: 


In, j , k] v - ^ ^ ^ X ‘ 1 y k 1 c n w (cos H) dx dy 


( 66 ) 


Substituting Equation 64 into Equation 61, 


C„ w ( cos H) = £ ( - 1 )« A" m (x - £)" C- m " 

m= 1 


|(y - J)” c;- i/2 (cos0) 


(67) 
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Substituting Equation 67 into Equation 66 and taking Equation 46 into account, 


n 

[n, j , k) = Y" (-l) m A y N l ' + m _ N* +m c 1 '' 1 ' 2 (cos^) 

v f v ' n ,m -j,n m,m k, n~m ,m m v ^ ' 


From Equation 53 we conclude that j and k must be integers of the form 


j - 2p - n , k - 2q-n, p,q - 0, 1, •" ,n 


( 68 ) 


Consequently 


C n v (cosH) 


y\ (n, 2p-n, 2q-n] v x 2 P* n y 2 ”’" 
p. q =0 


or 


C v (cos H) 

n v 


I> 


2p - n, 2q'- n] v cos [( 2p - n) 6 + (2q-n)0] 


(69) 


This last relation together with Equation 68 supplies the form of the addition theorem for 
Gegenbauer polynomials to be used in the expansions in the lunar theory. The polynomials 
[n, j , k] v are expanded in terms of Gegenbauer polynomials C m v ~ 1/2 (cos <£) . There is no actual ne- 
cessity to rearrange the right side of Equation 68 to represent an expansion in powers of cos<£. It 
is enough to provide a subroutine for the numerical evaluation of the polynomials C m v_1/2 (cosc£) 
for the given value of cos </>. 

The case 


v ~ \ » C n 1/2 ( cos H) - P n (cosH) 

requires a special consideration. In this case we have, from Equation 65, 


j , , k-+ n 

P n (cosH) = [n, j , k] 1/2 x j y k 

j , k-~n 


(70) 
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And by substituting 


P n m (cos 0) = (-l) m 2 m ^J, sin m ®C n ra _ + m I/2 (cos 0) 


= i"(i .) (*-J)"c;.r 




and 


P n m (cos </>) - ( 


-l) m 2 m m) sin m 0C n °L + m 1/2 (cos vt) 




,y + 71 

2 — 


into the addition formula 


P n (cosH) - P n (cos 6) P n (cos 0) + 2 ^ + ”y j p n m ( cos P n m ( cos cos m< £ 


we obtain 


P n (cosH) = £ A n i/2 ( x - ~) W C n ^ (v - y) m ^rT-V 72 (^T7 


(71) 


where 


A l /2 

n ,m 


- «-»■ (i »)’ 


/ 0 


A V2 
A n, 0 


Combining the result of Equation 71 with Equation 66, with v - 1/2, 


C n 1/2 (cosH) = P n (cos H) 
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I 


we obtain 


[n, j, k] U2 


E 

m — n 


hi/2W m * t/ J 2 N m + 1 '' 2 

n,m — j , n m , m ~k , n — m , m 


COS m 


(72) 


or, analogous to Equation 68, 


[n, j , k] , /2 


E 


A 1/2 N "i + l/2 j^m +1/2 
n , m “j , n~m, m - k, n~m , m 


m= 0 

where T m designates a Chebichef polynomial of the first kind, 
we conclude that 


j - 2p - n , k - 2q - n , p, q 


T m (cos <£) , (72') 

As before in the case of Equation 68, 

= 0, 1, ,n ; 


and consequently 


P n (cosH) 



2p - n , 2q ~ n] J/2 cos [( 2p - n ) 6 + ( 2q - n ) i/] 


(73) 


P. q-0 

Making use of the form of the addition theorem given in this section, we obtain the expansion 


m n 

(l + a 2 - 2a cos h) V = ^ * a n C ” (cos H) = ^ * a n ^ [n, 2p - n, 2q - n] v cos [( 2p - n) 6 + (2q - n)\//]. (74) 


This formula can be used to obtain the literal expansion of the disturbing function and its deriva- 
tives for the case of satellites in orbits highly inclined toward the orbital plane of the Sun. The 
polynomials [n, j, k] v can be expressed in terms of Tisserand polynomials Q[™) , which are ex- 
panded by Tisserand (1960) in powers of sin 2 1/2 and cos 2 1/2 . in the present work, the poly- 
nomials In, j, k] v axe expanded in powers of cos I, since this expansion is more convenient for the 
analytical treatment of the perturbations also in Delaunay 1 s canonical variables. 

If we decide to develop a purely numerical theory then we can also proceed in a different way: 
we can substitute the values e 0 , e 0 , Po > Po # > an< * 1 f rom the outset and use the formulas 


cos H - cos 2 2 cos (wj “ v' + w 4 j + sin 2 2 s in (w 1 + v ' - w 4 ^ 
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or the recurrence formula 


c ; 


(cos H) 


ti/2] 

- L 


c-iy 

m ! 


(v, j -m) 


(3 - m> ! 


' (2 cos H) 


j - 2m 


C/+ 2 (cos H) 


+ 2j + 2 


2v + j 

j + 2 


C. v (cosH) 


DEVELOPMENT OF THE DISTURBING FUNCTION 

In the theory presented in this work, the disturbing function has the form 

O - ^ v 2 (x 2 + z 2 ) tan 2 I - v 2 yz tan I + Q , (75) 

where, with properly chosen units, 



A = IV-V'I . 

In the system (x, y, z) the coordinates of the satellite are 

x = pcosWj , y = p sinWj , z - p£ , (77) 

and the coordinates of the planet in the same system are 

x' = p' cos (v' - w 4 ) , y ' = p' sin(v' -wj cos I , z' = ~p‘ s in (v ' ~ w 4 ) s in I . 

We obtain 

r • r # = pp' Jcos H - £ sin (v' - w 4 ) sin ij , (78) 

where 

cos H = cos Wj cos (v' - W 4 ) + sin Wj sin (v' - W 4 ) cos I . (79) 

We also have 

A 2 = (p 2 + p' 2 - 2p p' cos H) + |2pp' i sin (v' - w 4 ) sin I + p 2 £ 2 J ; (80) 
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and setting 


V = 1 + fi - cos H - 


= 2 £r i sin (v' - W 4 ) sin I + £ 2 

P 


we obtain from Equation 80 


i- 1 / 1 1 S 3 8 2 _ 5 8 3 

A ‘ y \ A o ■ 2 A 0 3 + 8 A 0 » " 16 A 0 7 


Substituting 


X 

v k+1 


2 

1 + ~p T C 1 k + 1/2 (cosH) +^C 2 h1/2 (cosH) + ••• 


into Equation 83 gives 


1 

A 


f-K + — C 1/2 + — C 1/2 ' + c 1/2 + c ,/2 ^ 
^ p' 2 C 1 p'3 C2 p'4 C 3 p's S J 


1- + c 3/2 + -El C 3/2 + -Bl. c 3/2 + ... 


2 \P' p' 2 


3 

4 w 3 


) 


+ l 82 + ti c i s/2 + 7T C / /2 + •■) 

-^ 83 (^ + TT c 1 7/2 + •••) + ••• • 


Taking Equation 78 in the form 


- P P 


c/ /2 - 2 P a 


8 + 


bp 2 t 2 


(81) 

(82) 


(83) 


(84) 


( 85 ) 


and substituting Equations 84 and 85 into Equations 76 and 75 gives 



where 


and 



0 ) 


i 


J. 

4 


v 2 p 2 ( 1 + cos 2w j ) tan 2 I 


- v 2 p 2 £ sin w 2 tan I 
+ -ip 2 £ 2 ^ 2 tan 2 I ' 


( 87 ) 


( 88 ) 


The quality $ 0 constitutes the main part of the disturbing function, and it is of the second order 
relative to the ratio of mean motions of the Sun and of the satellite; <t>j is of the second order 
relative to <P 0 . 

Let us introduce the set of the following auxiliary quantities and series which appear in the 
expansion of the disturbing function and in the expansion of its derivatives: 
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(93) 


r = r +p +p + 
1 1 1 11 1 12 1 13 


r i = r i0 + r ii + r i2 + ••• - 1 

> 1 , 

(94) 

8 0 = 2r 02 + 3r 03 + 4r 04 + •" ■ 


(95) 

8 1 = r n + 2r i2 + sr i3 + • ' • - 


(96) 

8 i = r il + 2P i2 + 3P i3 + •" ’ 

ill. 

(97) 

p p° r / t \ p 

= p" • ^[ 2sin ( v -w 4 )sinl • 

H • 

(98) 

A = M 2 (-5 r . + l < r 2-H e2r 3 + - 

-) • 

(99) 

* = P 2 ( 8 0 - \ e8 X + ! e28 2 ' B £30 3 

— ) • 

(100) 


These two last series are the basic ones in the theory presented here. The quantity p is analogous 
to the ratio of mean mean motions of the theory of Delaunay, and a is the parallactic factor of our 
theory. We obtain, after some easy transformations, 


p 0 °o = M 2 (r 0 er, + | € 2 r 2 - ^ € 3 r 3 + ---) . 


( 101 ) 


From the last equation and from Equations 91 through 98, we deduce 


'0 _ £_ 
5 o 


( et #„ ■ 


P d P P 0 

from which, taking Equations 19 and 21 into account, we obtain 


(102) 


(5$ r- ^ -i 

2 -j-1 = (l~g) 2 (l + s) (l + cos 2wj) tan 2 I - 2J sin Wj tan I + J 2 tan 2 I j 


- Si 2 


{* ■ £)' 
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Taking the two last equations and Equations 19 and 21 into account, we obtain 


p2™ = _£ 

p dp P„ 


_(**£• 7 


+ (1 - g) 2 ( 1 + s) (l + cos 2wj) tan 2 I - 25 sin Wj tan I + ^ 2 tan 2 i] 


M 2 t : 




’ (103) 


Now making use of 


dc; (x) 

—35— = 2pC n lV(*> 


we obtain, from Equations 90 through 94, 


ar\ 


d cos H 


p K o 

(2i + D- ' 7 ar ul 


i = 0, 1, 2, 


(104) 


Differentiating Equation 101 relative to w and taking Equations 79, 99, and 104 into account, we have 


*(Pq*o) 

dw 


2 [sin 2 sin (w x + v' - w 4 } + cos 2 -j sin ^ - v' + w 4 )] — aA 

Pn P 


(105) 


and from Equations 17, 21, and 88 

*K*,) 


<9w 


“ 2 ^ 1 - 4 ^ (1 + s) 2 (sin 2wj tan 2 I + 2^ cos tan i) . 


(106) 


From these two last equations we deduce 


/Q 2 a<t> 
p 0 <9w 


= { po ) { 4A [ sin2 \ sin K + v' -w 4 ) +COS 2 \ sin(w 1 -v‘ +w 4 )] ~ ■ ^ a 


^ ^ ( 1 + s ) 2 (sin 2wj tan 2 I + 2 £ cos w t tan *»}■ (107) 
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Substituting 


d (po^o) 


~ar 


= A , 


de 

TZ 


p * o 

Po ’ p' 


2 sin (v' -w 4 ) sin 1 + 2 


»o f' J 


into 


<?(p 0 ®o) _ 5 (Po®o) de 

5 1 " de dt 


we have 


*(p 0 ®o) p po' r , , v . T _ p Po' " 

71“ = To ' ~ [ 2 sin(v - w <) sin 1 + 2 JT 0 ■ T 7 V 


From Equation 88 we have 

d (Po ®i) 

di 


Po 


£-v 2 p 2 sin Wj tan I + p 2 S^v 2 tan 2 I - ^ 7 ^ ^ . 


and, taking Equation 22 into consideration, 


*(p 0 *i) 


as 


^ 1 ~ - ~ iP (1 + s) 2 (-sinWj tan I + £ tan 2 I ) - p 2 



Combining the last equation with Equation 108 we obtain 


d (p 0 *) 


as 


= 2 


p£ • ^[sin(v'-w 4 )sinl + 


+ ^ | j. ^ (1 + s) 2 (-sinWj tan I + £ tan 2 I ) ~ p 2 


and finally 


H = 2 (^) T 7 A [si n (v'- w 4 )si nI • ^7 t] 


+ (l~g) 2 (l + s) (-sin Wj tan I + £ tan 2 i) ~ p 2 
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The limit of the lower index of Gegenbauer polynomials in Equations 90 and 91 depends 
upon the values of a and £. If there is a necessity for developments for higher indices, we can use 
Equation 69 for the trigonometrical expansion of Gegenbauer polynomials in terms of the gener- 
alized Cauchy numbers. 

INTEGRATION OF THE DIFFERENTIAL EQUATIONS 

As the final output we must obtain the expansions of Brendel f s coordinates s, \p 9 and £ in Fourier 
series in four basic arguments: w n w 2 , w 3 , w 4 . The arguments w 1? w 2 , w 4 have a direct geometric 
meaning in the "nearly ideal" system of coordinates, which we have discussed previously. How- 
ever, w 3 requires additional discussion. The argument w 3 , which enters into the expansion of s, 

0, and £, is intimately connected with the argument v'. The angle v' defines the position of the 
Sim in its ellipse, and it is this "elliptic" v' which enters into the expansions given in Equations 87 
through 103 of the disturbing function and of its derivatives. Thus, these developments contain 
five arguments: w 2 , w 2 , w 3 , w 4 , and v'. We can set 


v' = w 3 + Sv' , (42 T ) 

where again Sv' is a series in four arguments, and we can eliminate v' in favor of w 3 by expanding 
0> o and its derivatives into a Taylor series relative to Sv'. However, in this expansion we must 
distinguish between w 3 as an argument in the series representing Brendel T s coordinates on 
the one hand, the "elliptic" values of v' and w 3 in Equation 42 T on the other. We have to eliminate the 
"elliptic" v' and, because we have two types of w 3 , we must do so by making use of two operators. 
One is the replacement operator R which means the formal replacement of the "elliptic" v' by w 3 : 

RF(w r w 2 , w 3 , w 3 ; v') = F(w, , w 2 , w 3 , w 4 ; w 3 ) . 

The other is the Taylor operator 


T = 



which transforms the function f(w 19 w 2 , w 3 , w 4 , v') into the function f(w 1? w 2 , w 3 , w 4 , v' + Sv'): 

TFjwj, w 2 , w 3 , w 4 , v') - F^Wj, w 2 , w 3 , w 4 , v' + Sv') 

Making use of the "product" 


K - RT 
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of both operators, we have from the last equation: 


KF^Wj, w 2 , w 3 , w 4 , v') - RF(w 1? w 2 , w 3 , w 4 , v' + Sv ' ) 

= F(wj, w 2 , w 3 , w 4 , w 3 +8v') = F(wj, w 2 , w 3> w 4 , v') 

+ 00 + 00 , V 

V 1 Sv' k <j k F _ Sv' k d k F\ 

*L k! <?v ,k L k! R \dv ,k / 


The operator R is analogous to the "bar-operator" of Hansen's lunar theory. Thus 


V" 1 5v' k d k 

Z_. k! R dv' 


k = 0 


(109) 


and Equations 24, 30 T , and 39" are rewritten in the form 


^ - a + b-pk(p^- 


d 2 s 
dv 


(30") 


d 2 j 
dv 2 


+ i + 2g(l - g) sin w 


= p + T K [( 1 + 5 2 )H- S H?)] 


(39") 


&P _ 2p 2 <9<J> 


dv 


Kav + d-g) 


[s ' 4(1 -x/i) ^]sinw tan I 




(24") 


where the expansion (109) of K is to be used. The expansion of Sv' in terms of the basic arguments 
and the value of m are determined from Equations 42 and 44. 

In Equations 30", 39", and 24", all arguments are now linear in v, and this facilitates the 
problem of integration. As the method of solving these equations, we suggest the process of 
iteration, partly because it makes the program homogenous and partly to avoid the question of 
convergence relative to the parameters e Q , e 0 ' , a, m, and I. We can substitute their values from 
the outset, before we start the process of iteration. At each iteration step, the frequency g of the 
argument w t is so determined that no term of the form A sin w x is present in £. In a similar 
manner we determine the frequency c of w 2 , by requiring that the only term of the form A cos w 2 in 
s be e 0 cos w 2 . The frequency of w 4 is simply (1 - g)/cos I . 
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EXPANSION OF THE RECTANGULAR COORDINATES 

The theory of satellites in Brendel's coordinates also favors expansion of the rectangular 
coordinates in the inertial system (xyz) into Fourier series. We have 


but 


r - p cos w J +Qsinw 1 +^Rj , 

P 0 ( 1 “ 0) 

P ~~ 1 + s 


[P, Q, R] 


COS W 4 

s in w 4 
0 


-cos I sin w 4 
cos I cos w 4 
sin I 


sin I s in w 4 
-s in I cos w 4 
cos I 


and we have finally 


x - p [_cos 2 ~2 cos |wj + w 4 ) + sin 2 cos ( w j ~ w 4 ) + C sin I siaw 4 

y - p [cos 2 j sin ^ + w 4 ) - sin 2 ^ sin (Wj - w 4 ^ - ^ sin I cos w 4 
z = p(sin Wj sin I + ^ cos I ) 

Substituting here the expansion of p, we have the expansions of the form 


(HO) 

(HI) 

( 112 ) 


X 


E< 


l l l 2 l 3 *4 


(*. 


Vj + i 2 


v + i v 
2 x 3 


3 + *4 V 4) 


(113) 


y. z 


L 


s. 


sin (i 1 v 1 + i 2 v 2 + i 3 v 3 + i 4 v 4 ) 


(114) 


where we set, as suggested by Equations 110 and 111, 


+ V, 


In these last expansions, three arguments v lf v 2 , v 4 , have nearly equal frequencies. 
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COLLECTION OF FORMULAS 


S - cos 2 ~2 cos ( w ! - v ' + w 4 ] + sin 2 cos ( w i + v ' - w 4 ) 
S 2 , S 2 , ••• 

F = i 

o - * 

P 

Q = FG 
GQ = Q x 
GQ 2 = Q 2 

GQ3 = Q 3 


L - 2 sin (v ' - w 4 ) sin I + Q 

e n £ LQ 

M = e + a£ 2 Q 

T = (l“g) 2 (l + s) (-sin Wj tan I + £ tan 2 i) 

U = sin(v' “ w 4 ) sinl + £Q 

W - ~(l~g) 2 (l"0)(sin2w 1 tan 2 I + 2£sinw 1 tanI^ ( a 2 £Q 3 

1 + s d 2 4 1 3 1 ds d <jj_ 

A 1 - 4> ^2 *2 1 ” 0 dv dv 
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B = 


3 15 M 3 1 + s td. p\ 2 

-J t 2 + ~g- 1 ‘ ’ - 2 (1 -^2 Vdv/ 


+ (l-g)(l+s)(2^-^ prj 37 j cos w tan I 


d£ 1 i 


) 


1 1 d£ . , (dl* 1 1 + £ 2 d/A 

2 dt dt + (1 -B) ^“dT " 2 1 - </> 37 j cos w tan I 


C = (1-g) 


[s 37 " 4 (! - 0) 37] sin w tan I 


^ 1 

C 1/2 = ^ 02 _ JL 

2 b 2 


S ^ 

r 1/2 = 4 ?2 - £ o 

S 2 b 2 s 


35 . 15 2 3 

C 4 1/2 = -g- S< - -4 S 2 + 8 


p 3/2 
'“'l 

= 3 S 

r 3/2 = 

^2 


r 3/2 - 

¥ S3 - ¥ s 

c 5/2 

= 5 S 

r 5/2 - 

^2 

35 g 2 __ 1 
2 b 2 

p 7/2 

= 7 s 

r 02 = 

= Q 2 C 2 1/2 
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r 03 = -Qs C 3 1/2 

r 04 = a2 Q4 C 4 /2 


®0 = ^02 + 3r 03 + 4P 04 + 

•"ll = Q! C 1 3/2 

r ,2 = a Q 2 C / /2 
r !3 = “ 2 Q 3 C 3 3/2 


p = p +p + r + 

1 1 * 11 12 13 


®X = r ,l + 2F 12 + 3F 13 + 


r 20 = G 


r 2 i = a Qj cp 2 

r 22 = « 2 Q 2 c 2 s/2 

P 23 = a3 Q 3 C 3 S/2 


p = p +p + r + r 

2 20 21 22 23 


®2 = r 21 + 2F 22 + 3r 23 + 


r 30 = G 
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r 3 , = a Ql c/^ 

r 3 2 = “ 2 Q 2 C/ /2 

r 33 = a3 q 3 c/' 2 


r 3 = r 30 + r 31 + r 3 2 + 


®3 ^3 1 + 2 ^3 2 + *^33 + 


a = m 2 {-- 2 r 1 + | £ r 2 --g e 3r 3 + ...) 


( 0 O - i «©1 + -| « 2 ® 2 - ll « 3 ©3 + '••) 


P 2 dp = f(MA + V) 


P 2 <?$ „ 

2 ?7^ = 4F 2 SaQA + W 


U'V 

P Jl = 2FQAU + T 

(l +Sv* -^7 t X Sv ' 2 — ~~~r + ■■■) 
\ dv 2 dv’ 2 ) 

.s = £ („J?)] 


d 2 i 

^ 2 + £ 2g(l-g)sinw - P 


ri + c 2 

/ 

£ ( 2 d<s > ) 

[_ 1 + s 

\ p di) 

1 + s V° dp )_ 


*L = _ 2 p1„™. c 

^ p 0 K dv + { ~ 


const, part in y p j/l + — - (l-i/>) 3 
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CONCLUSION 


Sv ' 


^period, part 



( 1 - 40 - 3/2 



^ P 0 3/2 ( 1 - 0) 3/2 (1 + s)" 2 


We have suggested a numerical theory of satellites which gives the expansion of the rectangular 
coordinates in terms of four arguments, using the expansion of Brendel’s coordinates as an inter- 
mediary step. This theory, like Hansen’s, can be used for highly inclined orbits up to the critical 
inclination, and it makes use of a "nearly ideal" system of coordinates. Unlike Hansen, we do not 
make use of the rotating ellipse as a reference orbit, and we have discarded the use of the w-function. 
Thus, the theory given here leads to the numerical expansion of the rectangular coordinates in a 
more direct way. Recently, the author has suggested a modification of Hansen’s lunar theory which 
was programmed and compared with Hansen’s theory to the extent and to the accuracy necessary 
for computing ephemerides of the satellites of outer planets (1 x 10 -5 in r, l°x 10” 4 in angles). Our 
next goal is the programming of the theory presented here and comparison of the expansion with 
those given by other theories. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, January 10, 1967 
188-43-01-01-51 
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Appendix A 

Symbol List 

the mass of the Sun, the mass of the planet, and the gravitational con- 
stant are set to be one. 

the inertial system of coordinates, with the X- and Y- axes in orbital 
plane of the Sun. 

the mean inclination of the orbital plane of the satellite toward the 
orbital plane of the Sun. 

the unit vector directed from the planet toward the ascending node of 
the mean orbital plane of the satellite. 

the unit vector normal to the mean orbital plane of the satellite. 

the position vector of the satellite with respect to the planet. 

the projection of r on the mean orbital plane of the satellite, 
the position vector of the Sun with respect to the planet. 

the semi-latus rectum of the solar orbit, 
the true longitude of the Sun. 
the projection of r on P. 
the projection of r on Q. 
the projection of r on R. 

the angular speed of rotation of the mean orbital plane of the satellite 
around the Z-axis. 

the projection of y on Q. 

the projection of 7 on 1 
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€y 


v 


W - Wj - gV+Cj 
W 2 = CV + C 2 

W 3 = mV + C 3 

w 4 = ( 1 - g) V sec I + C 4 


P 

Po 

-p 0 'P 


s 


the "nearly ideal" system of coordinates located in the mean orbital 
plane of the satellite and rotating with the angular velocity -vR rela- 
tive to the frame (P, Q). 

the polar angle of p in the system (€y). 

the mean argument of the latitude of the satellite. 

the mean true anomaly of the satellite. 

the mean longitude (relative to v) of the Sun. 

the longitude of the mean ascending node. 

the areal velocity of p in the system (£y). 

the mean areal velocity of p in the system (s y). 

the perturbations in the mean areal velocity of p relative to (£y) 

P = P 0 (1 

p 

the Brendel’s coordinate defined by the equation: p - 



<J> = (x 2 + z 2 ) - \vyz + H 


the disturbing function associated with the solar perturbations in the 
motion of the satellite. 

the modified disturbing function. 

the Cauchy number. 

the generalized Cauchy number. 
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